Accessing the Neel phase of ultracold fermionic atoms in a simple-cubic optical lattice 
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We examine the phase diagram of a simple-cubic optical lattice half-filled by two species of fermionic atoms 
with a repulsive s-wave contact interaction. We use the Hartree approximation in the regime of weaker interac- 
tions, and a Hubbard model approximation in the Mott insulating phase. The regime where the Neel phase of 
this system is likely to be most accessible to experiments is at intermediate lattice strengths and interactions, and 
our two approximations agree fairly well in this regime. We discuss the various issues that may determine where 
in this phase diagram the Neel phase is first produced and detected experimentally, and analyze the quantum 
phases that occur in the intermediate lattice strength regime. 
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One of the next notable milestones in the production of new 
strongly-correlated states of ultracold atoms is expected to be 
the antiferromagnetic Neel phase of two hyperfine species of 
fermionic atoms in an optical lattice . Important progress 
towards this goal includes the recent realization of the Mott 
insulating phase with fermions IHQ], and the demonstration 
of controllable superexchange interactions in an optical lat- 
tice, albeit with bosonic atoms [5]. The interest in experimen- 
tally realizing this Neel phase is partly because it is a natu- 
ral next step in developing an ultracold-atom quantum "sim- 
ulator" of the fermionic Hubbard model with repulsive inter- 
actions, which is one of the prototypical models in strongly- 
correlated quantum condensed-matter physics. As we will ar- 
gue, if the objective is to access a phase with antiferromag- 
netic order, then it may be favorable to explore the regime of 
an optical lattice of intermediate depth, where the effective 
Hamiltonian is more than just the standard one-band Hubbard 
model. The low temperature phase diagram in this regime 
is worth exploring in its own right, as a new, interesting, and 
hopefully experimentally accessible quantum many -body sys- 
tem. 



An outline of our paper is as follows: first we review the 
physics of two-species fermions in an optical lattice, with 
the interactions tuned by a nearby Feshbach resonance. We 
then describe the Hartree approximation which captures the 
ground-state phase diagram, and use this calculation to esti- 
mate the line in the phase diagram along which antifeiTomag- 
netic interactions are strongest. For another approach to parts 
of the same phase diagram, we map the system onto a lattice 
Hubbard model, which is appropriate for a strong lattice po- 
tential, and study which terms dominate the Hamiltonian as 
one weakens the lattice away from the strong lattice regime. 
We also estimate the Neel temperature in this strong-lattice 
approach, and show that it is maximized along a line in the 
phase diagram that agrees well with the line of maximal ex- 
change interactions estimated using the Hartree approxima- 
tion. Finally, we discuss experimental issues such as equi- 
libration, Feshbach molecule formation, and the size of the 
Mott phase in a trap. 



THE MODEL 

We consider the model that experimentalists are focus- 
ing on, which consists of two hyperfine species of fermionic 
atoms of mass m in a simple-cubic optical lattice with lattice 
spacing d and a single-atom optical potential of the standard 
form [il: 

+ sin — 
d 



Vi{x, y, z) = Vo{sm^ — + sm^ -f + sm^ — ) 



(1) 



Atoms of opposite species interact repulsively with a contact 
(s-wave) interaction, whose strength can be adjusted by tun- 
ing the applied magnetic field relative to a nearby Feshbach 
resonance, and can be quantified by the atom-atom scattering 
length a.; > in the absence of the lattice. We call these two 
states "up" and "down" and treat them as the two states of 
a spin- 1/2 degree of freedom. The s-wave repulsive interac- 
tion is only between atoms of opposite spin. To lowest order 
in (fls/d), this 2-atom interaction is the standard regularized 
contact potential [12], 

^2f|-rx) = d(ri-ri)—r, (2) 

m or 

where r is the distance between the two atoms. We will work 
in the thermodynamic limit of this model at half filling, and 
we are particularly interested in the Mott insulating phase. 

This system is appealing because it behaves effectively like 
a one -band Hubbard model at large Vq, with on-site repulsive 
interaction U and nearest-neighbor hopping matrix element 
—t, so that experiments may eventually resolve long-standing 
questions about the Hubbard model. However, as Vq is low- 
ered one eventually reaches an intermediate lattice-strength 
regime where higher bands can no longer be neglected. To 
address this portion of the phase diagram, we use the self- 
consistent Hartree approximation, which includes admixtures 
of the single-atom states in higher bands. 



HARTREE APPROXIMATION 

When considering the regularized contact potential in 
mean-field theory, there is no exchange term, implying that 
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the Hartree and Hartree-Fock approximations are identical 
here. The total effective potential "seen" by the atoms with 
Sz — —<J in the Hartree approximation is thus 

V[^ff\r) = + V,{r) , (3) 

m 

where (r) is the number density of atoms with Sz = cr at 
position r. 

For each point in the ground-state phase diagram, speci- 
fied by the lattice intensity Vq and the interaction as, we solve 
the Hartree equations numerically by discretizing them in mo- 
mentum space, and iteratively achieving self-consistency. We 
obtain up to 3 different low energy self-consistent Hartree 
many-body states at density one atom per lattice site, and de- 
termine which state is of the lowest energy and thus is the 
Hartree ground state. In the paramagnetic state we impose the 
constraints (f) — ni (r) = (r + S) at all r, where S is any 
nearest neighbor vector of the simple cubic lattice. This para- 
magnetic Hartree state exists for all values of as and Vq, since 
this constraint is preserved by the iterations. In this state the 
lowest Hartree bands are always partially occupied, so there 
are Fermi surfaces. For the Neel antiferromajnetic state we 
impose only the constraints n|(r) = n|(r + (5), thus permit- 
ting two-sublattice Neel order At weak enough lattice and 
interactions, the only self-consistent Hartree state is in fact a 
paramagnetic state, i.e. it has no magnetic order. When an 
ordered Neel state exists, it can be either Mott insulating, with 
the lowest Hartree bands full and a Mott-Hubbard gap to the 
next bands, or the two lowest Hartree bands for each species 
can overlap in energy and each be partially occupied, with 
Fermi surfaces. We also look at ferromagnetic states that have 
different densities of f and j atoms but do not break the dis- 
crete translational symmetries of the lattice. Again, there is a 
portion of the phase diagram where there is no self-consistent 
ferromagnetic state. And when there is such a state, it can 
be either a band insulator or have Fermi surface(s), and it can 
be either fully or partially spin-polarized. In the latter case 
(PPF in Fig. 1) there are 3 partially occupied Hartree bands. 
In these candidate Hartree ground states, the fermions occupy 
only the lowest Hartree bands, but in terms of free noninteract- 
ing fermions, we have included states extending out to many 
Brillouin zones and thus the Hartree states are admixtures of 
multiple bands of the noninteracting system. Specifically, for 
the parameters in Fig. 1, a 20 x 20 x 20 grid of momentum 
points in each Brillouin zone and a 9 x 9 x 9 grid of zones 
was enough to achieve convergence everywhere. 

The resulting ground-state Hartree phase diagram of this 
system is presented in Fig. 1. At any lattice strength, the 
paramagnetic Fermi liquid phase exists at weak enough in- 
teraction. At strong interaction, the Hartree approximation al- 
ways produces a ferromagnetic ground state, due to the classic 
Stoner instability. The ferromagnetic phase is a band insula- 
tor for Vq > 2.2Er, with a band gap above the filled lowest 
fully spin-polarized band. For Vq < 2.2Er the lowest two 
bands for the majority-spin atoms overlap and there is instead 
a ferromagnetic Fermi liquid. We show this ferromagnetic 




FIG. 1 : Ground-state phase diagram for filling one fermion per lat- 
tice site in the Hartree approximation. The phases shown are the 
antiferromagnetic Mott insulator (AFI); paramagnetic (P), antiferro- 
magnetic (AFM) and partly- (PPF) and fully-polarized (FM) ferro- 
magnetic Fermi liquids; and the ferromagnetic band insulator (FI). 
The solid line marked J™"^ indicates where the Hartree estimate 
of the effective exchange interaction J as a function of the lattice 
strength Vo is maximized for each interaction as , and the dashed line 
near it shows where our estimate of the Neel temperature Tjv gets 
maximized under the same prescription (see text). The dash-dotted 
U = 14t line is near where the DMFT estimate of the entropy is 
maximized at the Neel ordering temperature Tjv [§]. The square 
on each line denotes the overall maximum of J or Tn along that 
line. The lattice intensity Vo is given in units of the recoil energy 
Er = {■Khf/{2m(f). 



part of the phase diagram for completeness, but it is impor- 
tant to emphasize both that the Hartree approximation is not 
to be trusted at such strong interactions, and that systems of 
ultracold fermionic atoms are likely to be highly unstable to 
Feshbach molecule formation when brought this close to the 
Feshbach resonance. Thus we do not expect such a ferromag- 
netic phase to be experimentally accessible for these systems, 
even if it does exist in a model that ignores the instability to- 
wards Feshbach molecules. Finally, for lattices above a cer- 
tain minimum strength there is a Neel-ordered ground state 
at intermediate values of the repulsive interaction. This Neel 
state is a Mott insulator over most of the phase diagram, and 
becomes an antiferromagnetic Fermi liquid over a small sliver 
of the phase diagram at weak lattice and weak interaction be- 
tween the paramagnetic and Mott insulating phases. It is in 
this lower (smaller as) portion of the phase diagram where 
we believe the Hartree approximation is qualitatively correct, 
producing the paramagnetic and antiferromagnetic Fermi liq- 
uid phases and the Mott insulating Neel phase. 

To quantify the energy scale associated with Neel order- 
ing in the antiferromagnetic phase, we (crudely) estimate the 
nearest-neighbor antiferromagnetic exchange J by taking J/2 
to be the energy difference per bond between the Hartree anti- 
ferromagnet and ferromagnet (at low as and/or Vq, the ferro- 
magnet has zero magnetization, so is really the paramagnet). 
For each given interaction a^, we locate the lattice intensity 
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FIG. 2: Plot of the Hartree estimate of the antiferromagnetic ex- 
change coupling, J, as a function of Vb for Us/d = 0.08, compared 
with the estimate from the perturbative expansion at strong lattice. 
The red (full) line shows the Hartree estimate of J, while the blue 
(dotted) line gives the perturbative estimate, discussed in the text. 
They match very well at high lattice depth Vb . 



Vb where this Hartree estimate of J has its maximum, and 
indicate these maxima by the red (full) line marked J™"^ in 
Fig. 1. It is somewhere near this line that antiferromagnetic 
interactions are strongest, and the Neel phase survives to the 
highest temperature. The overall global maximum in this es- 
timate of J occurs at Qg/d ^ 0.15 and Vb = SE'r, where 
Er = {ttH)'^ / {2m(P) is the "recoil energy". This point is 
marked with a small square on the full red line in Fig. 1. 

As an example, in Fig. 2 we show our Hartree estimate of 
J vs. Vb at Us/d — 0.08, which is the largest interaction ex- 
plored experimentally in the study of ref. |3], although these 
experiments only reached temperatures far above those of the 
Neel phase. In Fig. 2 we also plot the estimate of J from a 
strong lattice expansion which we discuss in the next section, 
and it coincides well with the Hartree estimate for Vq > bEr- 

The line in Fig. 1 where the Hartree estimate of J is maxi- 
mal not only indicates roughly where Tn is maximized, but it 
also occurs at the location of a crossover in the nature of the 
antiferromagnetic state. In the regime of larger Vq and a^, the 
system is a "local moment" magnet: almost every well of the 
optical lattice is occupied by a single atom with a "spin", with 
very few wells either empty or multiply-occupied. Here we 
can make a Hartree ferromagnet as well as an antiferromag- 
net and compare their energies to obtain our estimate of J. In 
fact we could make many other metastable magnetic Hartree 
states with ordering at essentially any momentum in the Bril- 
louin zone. Here J is primarily a superexchange interaction 
J ~ t^/J7, and it increases as the lattice depth Vq is reduced 
(see Fig. 2), since this increases the hopping t between wells 
and reduces the interaction U between two unlike atoms in the 
same well. 

In contrast, the regime of the antiferromagnetic phase that is 
at lower Vq and Us near the paramagnet is not a local-moment 
regime, but instead is a spin density wave (SDW). Here one 
can make a magnetically-ordered Hartree state only at mo- 
menta near those that cause substantial Fermi-surface nest- 



ing, which is initially only near the corners of the first Bril- 
louin zone. In this regime we do not have a ferromagnetic 
Hartree state and our estimate of J is obtained by subtract- 
ing the energy of the Neel state from the paramagnet, so this 
J should not be interpreted simply as a spin-spin interaction. 
We now understand why this estimate of J increases with in- 
creasing Vb (see Fig. 2), since this increases the interaction 
that causes the magnetic ordering, and decreases the hopping 
that favors the paramagnetic Fermi liquid. The effective J is 
maximized at the crossover between the SDW and local mo- 
ment regimes. This crossover shows up in the metastable fer- 
romagnetic Hartree state as two closely-spaced phase transi- 
tions between un-, partially-, and fuIIy-poIarized. 



STRONG LATTICE EXPANSION 

We now proceed to mapping our system onto a lattice Hub- 
bard model, which is an appropriate description at deep lat- 
tice. This mapping uses the single-atom Wannier orbitals as 
the basis states uJ, Ig, uM- The standard one-band Hubbard 
model includes only the lowest-energy Wannier orbital at each 
lattice site as its Hilbert space, and neglects all terms in the 
Hamiltonian except for nearest-neighbor hopping and on-site 
interaction between different hyperfine states in the lowest 
band. 

The single-atom optical-lattice potential (1) that we con- 
sider is a separable potential, so the energy eigenstates of a 
single atom in this potential can be chosen to be the prod- 
uct of one-dimensional (ID) eigenstates along each direction. 
We solve for these ID bands, and use the gauge degree of 
freedom in the definition of Wannier states to minimize the 
spatial variance of their probability distributions. Thus we ob- 
tain the normalized wavefunctions u>„ (x) of the maximally- 
localized ID Wannier orbitals for each band n (n — 0,1, 2, ..) 
lfl3ll . The 3D Wannier orbitals are then if)n^nyn^ix,y, z) = 

WnAx)Wny{y)WnAz)- 

Since our focus here is on the Mott insulator with density 
one atom per lattice site, low temperature and fairly strong 
lattice Vq, the atoms will primarily be in the lowest band, 
rix — ny = Uz = 0. The nearest-neighbor hopping matrix el- 
ement Iq in this noninteracting single-atom band is a strongly 
decreasing function of Vq, the lattice strength, behaving as 

for large Vq fd] 



for large Vq [6]. The wavefunction 
wq (x) of the lowest Wannier orbital at a given lattice site has 
its maximum amplitude at a; = 0, the center of that well of 
the optical lattice, while its amplitude changes sign and is of 
smaller magnitude in the nearest-neighbor wells. 

The antiferromagnetic Mott insulator occurs in the regime 
where the on-site repulsive interaction U is stronger than the 
hopping t. The expectation value of this interaction energy for 
two atoms of opposite spin occupying the same lowest Wan- 
nier orbital is our first estimate of the strength of the on-site 
interaction UQUn^riii in the corresponding one-band Hubbard 
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model: 



(4) 



In the Hubbard model, when adjacent sites i and j are each 
singly-occupied by atoms with the same spin, then the hop- 
ping between those two sites is Pauli-blocked. When these ad- 
jacent sites are each singly-occupied by opposite spins, then 
virtual hopping between these sites, treated in second-order 
perturbation theory, allows them to lower their energy and 
thus generates the antiferromagnetic superexchange interac- 
tion JsiS, ■ Sj - I) with Js = U^/U. 

The leading corrections to the Hubbard model approxima- 
tion to our system in the regime of interest are due to the in- 
teractions between atoms of opposite spin occupying lowest 
Wannier orbitals on nearest-neighbor sites i and j. By ex- 
panding fully the interaction term in the Wannier basis, and 
evaluating all the terms, one can show that the most important 
term in limiting the Neel phase is the "direct" interaction ISJ 



dxwo(x)]^ / dywl{y)wl{y + d) (5) 



between atoms of opposite spin in adjacent orbitals. This term 
is due to the overlap of the probability distributions of adja- 
cent Wannier orbitals. It raises the energy of the Neel state, 
and thusjjroduces a direct ferromagnetic exchange interac- 
tion Jf{Si ■ Sj — i) with Jf — — 2[/„„ < that partially can- 
cels the antiferromagnetic superexchange Jg that occurs in the 
Hubbard model. It is primarily this ferromagnetic interaction 
that stops and reverses the increase in the net antiferromag- 
netic interaction J as one moves towards stronger interaction 
and a weaker lattice while staying near the values of the lattice 
strength Vq that maximize J. At the global maximum of J, 
indicated in Fig. [H we find J/ ~ — Js/4. Another (small) ef- 
fect of Unn that we include is that it changes the effective on- 
site Coulomb interaction in the antiferromagnetic Mott phase, 
leading to an effective interaction U ~ Uq — 6J7„„. 

For large enough this direct ferromagnetic exchange is 
stronger than the superexchange and thus we have the ground- 
state phase transition to the ferromagnetic phase, as indicated 
in Fig. 1. We see in Fig. 1 that the line where | J/| = Jg 
is close to the line where the Hartree calculation also signals 
this phase transition. If this transition is real, further-neighbor 
interactions might cause other magnetically-ordered phases to 
occur near it in the phase diagram, but this is an issue that 
we do not explore here, since we expect that this portion of 
the phase diagram will not be accessible to cold-atom experi- 
ments, as discussed above. 

The next leading correction to the Hubbard model is an ad- 
ditional hopping term of the same sign as to, sometimes called 
correlated hopping in the literature ||6[], 

ti = f dxw^ix)]^ [ dywl{y)wQ{y + d) (6) 



when the two sites are each singly-occupied by opposite spins. 
The effective hopping that enters in the superexchange pro- 
cess at this level of approximation is thus t — to + tj. 



More generally, the hopping matrix element for an atom with 
Sz = — cr hopping between nearest-neighbor sites i and j is 
t ^ tQ + n^tj, where here is the total number of atoms 
with 5^ = (T on sites i and j. 

Summarizing, once we include these leading effects due to 
the nearest-neighbor interaction, the effective Hamiltonian in 
the vicinity of the antiferromagnetic ground state of the half- 
filled Mott insulator has hopping t — to + tj, an effective 
on-site interaction U = Uo — 6?7„„, and an additional ferro- 
magnetic nearest-neighbor exchange interaction Jy = —2Unn 
when both sites are singly-occupied. [We note in passing that 
there are other condensed matter systems which behave effec- 
tively like a Hubbard model with ferromagnetic interactions, 
such as V2O3 ifTsll.l We can use this approximation to pro- 
duce an estimate of U / 1 for each point in the Mott insulating 
regime of our phase diagram, and we indicate in Fig. 1 the 
lines along which U /t = 6 and U /t ~ 14. 

Now that we have our effective model, we use it to deter- 
mine which values of the parameters will maximize the ro- 
bustness of the antiferromagnetic Mott phase. To that end, we 
start by looking at the Quantum Monte Carlo simulations lll411 
of the simple-cubic fermionic Hubbard model, which show 
that for a given U, the critical temperature T/v for the Neel 
phase is maximized near U/t ~ 6. From Fig. 1, we learn that 
the line along which U = 6t, as estimated from the Hubbard- 
model approximation is in reasonably good agreement with 
the location of the maximum J as estimated in the Hartree 
approximation. Thus we can hope to start from the QMC re- 
sult for T/v, and study perturbative corrections to it using the 
strong-lattice expansion, to obtain a reasonable approximation 
to Tn four our model. 

To estimate the Neel temperature of our system we propose 
the following approximation: For the Hubbard model with- 
out Jf, we have estimates of its Neel tem pera ture Tj^ it,U) 
from quantum Monte Carlo simulations |I14'1. In the regime 
of interest, this Neel ordering is due to the antiferromagnetic 
superexchange interaction Jg — Ai^ /U between neighboring 
singly-occupied sites. When we include J/ < the magnetic 
interaction is reduced, and we will approximate the resulting 
reduction of T/v as being simply in proportion to the reduction 
of the total nearest-neighbor exchange interaction: 



TN{Vo,as) 



it,U). 



(7) 



This approximation is certainly accurate for large U/t, where 
the system is well-approximated by the Heisenberg antifer- 
romagnet. By using this approximation at the maximum of 



( H 

Tj^ (t, U) near U = Qt, we are assuming that the reduction 
there by roughly a factor of 2 of TV from its Heisenberg value 
lITill is mostly due to the dilution of the Heisenberg antiferro- 
magnet by nonmagnetic empty and doubly-occupied sites and 
not due to a large change in the superexchange interaction be- 
tween singly-occupied sites. 

Using eq. (7), we can maximize Tm as a function of Vq 
for each value of a^/d. The resulting "optimal" values of Vq 
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T/v at the optimal V 

Tj^{Er) 
0.04 r 




FIG. 3: Solid line: our estimates of the maximal Neel temperature, 
Tm, as a function of as/d. For each value of as, Tn is maximized 
by varying the lattice depth Vo- Dashed line: we also plot Tn at the 
line U = 14t, which is near where the critical entropy is maximized 
in DMFT H. 



are indicated by the dashed red line in Fig. 1. The line of 
maximal Tn we obtain this way is quite close to the full red 
line along which the Hartree estimate of J is maximized. The 
global maximum of this estimate of Tn (marked by a small 
square on the line in Fig. 1) also coincides well with the global 
maximum of the Hartree J. Thus we see that in this regime 
these two approximations agree rather well, suggesting that 
these estimates are indeed to be trusted. We also show the line 
along which U — 14t; this is near where DMFT calculation 
show the critical entropy S{Tn ) to be maximal (@] . Along this 
U = 14i line the maximum in T^ also occurs near Qs/d = 
0.15. 

In Fig. 3 we show Tn as a function of as/d at the value of 
Vb that maximizes our estimate of Tn, as well as at the value 
of Vo that gives U = \At and thus is near the maximum of 
S{Tn) as estimated in DMFT |6]. Note that in Fig. 3 the 
horizontal scale for as/d is logarithmic, so Tn drops rather 
weakly as Og is reduced. This means that if the optimal value 
of as can not be reached due to Feshbach molecule formation, 
this may not "cost" very much in terms of the suppression of 
Tn- On the other hand, increasing U / 1 can produce a substan- 
tial decrease in Tn- it drops by about a factor of two between 
its maximum at [/ = 7t and U = lAt- For ^Li and d ^ 500 
nm, the overall maximum of Tn is a temperature of about 40 
nK. 

The strong lattice approximations we have used here are 
those appropriate for the antiferromagnetic Mott insulator, 
and are based on the inequalities on energy scales Eq > U > t, 
where Eq is the expectation value of the single-particle energy 
in a lowest Wannier orbital. By calculating higher order terms, 
we can check the reliability of the estimates coming from this 
strong lattice expansion. We do find that this expansion is be- 
ginning to break down in the vicinity of the parameter values 



that maximize J. Thus, although we expect these approxima- 
tions to give reasonably reliable rough estimates of the maxi- 
mal values of Tn, there are many higher-order effects that we 
are ignoring that may alter these estimates (our calculations 
suggest on the ~ 20% level). At the maximum of Tn, \Jf\ 
is about 25% of Js- The correction to Jg due to tj is also 
of roughly this size, but its dependence on as is much weaker, 
which is why Jf is the important actor in setting the maximum 
in Tn- 

We have analyzed in perturbation theory many corrections 
beyond those included above. We find that at the maximum 
of Tn (both the global maximum and the maximum along the 
U = 14t line) the strongest next coiTection is the fourth-order 
process illustrated in Fig. HI it alters the effective J by about 
10%. 




(1) 

FIG. 4: The strongest higher-order process contributing to the energy 
of the antiferromagnetic Mott insulator at the maxima of Tn shown 
in Fig. [T] It consists of (1) a nearest-neighbor hop in the lowest (S) 
band, (2) an on-site "pair hopping" of both atoms up to the next (P) 
band, (3) on-site pair hopping back to the S band, and (4) a nearest- 
neighbor hop back to the original configuration. At both maxima of 
Tn, this process corrects Js by about 10%. This indicates that mix- 
ing with higher bands is significant, which we don't capture with our 
strong lattice expansion, but do capture with our Hartree calculation. 



TABLE I: The values of the various energies at the two Tn maxima. 
The " Ja correction" corresponds to the process detailed in Fig. |4] 





Global maximum 

ofTjv 


Maximum 

of Tn SLtU ^ Ut 


eo (Er) 


4.2 


5.5 


Uo { Er ) 


0.9 


1.3 


to (Er) 


0.11 


0.07 


tl (Er) 


0.02 


0.02 


Js (Er) 


0.08 


0.03 


Jf (Er) 


-0.02 


-0.008 


Js correction (Er) 


0.007 


0.004 



We have seen that our two approaches give consistent re- 
sults in the region of intermediate lattice strength. We now 
briefly discuss possible improvements on what we report here. 
Since our perturbatively-based approximations may be break- 
ing down near this regime of interest where J is maximized, 
it would be nice to have a more systematic approach that can 
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obtain more precise and reliable estimates of the phase dia- 
gram in this regime. For example, quantum Monte Carlo sim- 
ulations might be possible at this maximal Tjv, although of 
course the famous fermionic "minus signs" may prevent this 
from being feasible in the near term. Another approach would 
be a systematic correction of the Hartree calculation. 

Since our calculations are in the thermodynamic limit, we 
are assuming that finite-size effects do not substantially shift 
the transition at TV or alter the ground-state phase diagram. 
The experiments use clouds of well over 10'' atoms [3], so this 
should be true unless a trap of extreme aspect ratio is used that 
reduces a linear dimension of the Mott domain to well under 
lOrf. 



EXPERIMENTAL CONSEQUENCES 

We are now in a position to discuss the experimental con- 
sequences of the calculations above. The assumption usu- 
ally found in the literature is that experiments will be able to 
evolve the system adiabatically, and that the natural variable 
is therefore the entropy, instead of the temperature. We would 
like to point out some caveats to this. 

The system must be cooled to low temperature T and thus 
low entropy S, and equilibrated at the point in the phase dia- 
gram that is being measured. This cooling may be done un- 
der some other conditions (e.g., with the lattice turned off 
and other parameters optimized for cooling), with the "pre- 
cooled" system then moved adiabatically to the conditions 
of interest for measurement Hi 0, S Isl For this to 
work, the time scales of the system must be such that this 
can be done without strongly violating adiabaticity. For the 
Neel phase of the Mott insulator, this will limit how small the 
hopping energy t and the antiferromagnetic superexchange in- 
teraction J can be, since the system must be able to remain 
near equilibrium, adiabatically rearranging the atoms so that 
there is one atom per lattice site and these atoms are antifer- 
romagnetically correlated. Alternatively, the system might be 
actively cooled under the conditions of measurement, but this 
again requires the system to be able to equilibrate under those 
conditions. Thus, either way, this constraint limits equilib- 
rium access to the strong-lattice portion of the phase diagram, 
where the exchange J and/or hopping t are too slow to allow 
equilibration and adiabaticity. Thus the Neel phase is going to 
be most accessible to experiment in some regime of interme- 
diate lattice strength Vq. These considerations may also limit 
how large the interaction can be made, since strong repulsion 
tts suppresses the superexchange rates and thus can limit spin 
equilibration in the Mott insulating phase. 

Another constraining issue is that as one approaches the 
Feshbach resonance from the repulsive side, two atoms will 
interact repulsively only as long as they scatter in a way that 
is orthogonal to the molecular bound states. Thus the atoms 
must remain metastable against forming Feshbach molecules. 
In the absence of the lattice, the rate of molecule formation 
grows as ^ as the Feshbach resonance is approached 1, 1 1.1 . 



This will limit how large the interaction ag can be made, in 
a way that is presumably less of a constraint as the lattice is 
made stronger so that it keeps the atoms apart. 

As we mentioned earlier, the DMFT estimate of the max- 
imum of S{Ti\f) is near U — 14t |6]. In fact, the maximum 
in S{Tn) is rather weak, being only a little higher than that 
of the Heisenberg limit U ^ t. Thus one might more use- 
fully say that the critical entropy is nearly maximized for any 
U > 12t. As we can infer from Fig. I, U > 12t occurs at 
relatively large Ug and Vq, and the issues raised above may 
force experiments away from the maximal entropy line, and 
towards the line where the exchange interactions and T/v are 
maximized. 

So far we have assumed an homogeneous system, but the 
trapping potential is actually nonuniform, which leads to inho- 
mogeneity of the local equilibrium state in the trap. The spa- 
tial size of the region occupied by the antiferromagnetic Mott 
phase will increase with increasing U /t, since U increases the 
Mott "charge" gap. This effect means larger U/t should favor 
detection of the Neel phase; of course the optimal U/t will be 
some compromise between this and the other issues discussed 
above. Using our Hartree calculation to obtain the charge gap 
in the Mott phase, we find within the local-density approxi- 
mation (LDA) that at the point in the phase diagram where J 
is maximized, the Mott phase should occupy about half of the 
linear size of the trap. This is encouraging, but the Hartree 
approximation likely overestimates the Mott gap, as the true 
gap should be renormalized downwards by spin fluctuations. 

Finally, we come to novel experimental predictions that re- 
sult from our calculations. The Hartree calculation predicts 
that there is interesting and as yet unexplored physics at low 
to intermediate lattice strength, and weak coupling. We men- 
tioned already that the effective model at around optimal val- 
ues of Vq and is that of a Hubbard model with ferromag- 
netic correlations, which is interesting in its own right. The 
different phase boundaries that we uncovered may also be 
worth exploring. For example, the quantum phase transitions 
between the Mott and metallic Neel phases and the paramag- 
net occur in parameter regimes that are quite accessible to the 
experiments, although it may not be possible to see their ef- 
fects at accessible temperatures, since T/v decreases strongly 
as this weak-coupling regime is approached. 

CONCLUSION 

We have shown that to maximize antiferromagnetic interac- 
tions for fermionic atoms in an optical lattice one must explore 
the regime of intermediate lattice depths, where the system 
has significant deviations from the standard one-band Hub- 
bard model. We have found that the nearest-neighbor direct 
ferromagnetic exchange is the most important correction to 
the Hubbard model that limits the maximal exchange ,/, and 
therefore the maximal Neel temperature ■ 

There are also higher-order corrections to the Hubbard 
model: virtual hopping into higher bands and other higher- 
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order processes. The relative contribution of the higher-order 
corrections in the vicinity of the optimal J drops exponen- 
tially as one goes to smaller interaction and thus a larger 
Vq. Thus our perturbative calculation should yield accurate 
results in the large Vb (strong lattice) regime. 

We included a subset of all higher-order corrections via a 
Hartree calculation, and used it to find an estimate of the line 
where the exchange J is maximized. This line coincides well 
with the line where our estimate of Tm is maximized, obtained 
by using quantum Monte Carlo results and the strong lattice 
expansion. 

For quantitatively more accurate results in the intermedi- 
ate lattice depth regime, one needs to resort to more system- 
atic quantum calculations. Of course this is a system of many 
fermions, so it is not clear whether this regime can be accu- 
rately treated in some form of quantum Monte Carlo simula- 
tions. 

Experiments on these systems still need to reduce the tem- 
perature by a substantial factor before they are able to access 
magnetically ordered phases of fermions in optical lattices. 
Once they bridge this gap, our results above suggest that the 
Neel phase will be most accessible in the intermediate lattice 
strength regime. If that is the case, then this "quantum simula- 
tor" should be able teach us about more than just the standard 
one-band Hubbard model. 
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